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Summary 



9 Ray-finned fishes constitute the dominant radiation of vertebrates with over 30,000 species. 

10 Although molecular phylogenetics has begun to disentangle major evolutionary relationships 

11 within this vast section of the Tree of Life, there is no widely available approach for effi- 

12 ciently collecting phylogenomic data within fishes, leaving much of the enormous potential 

13 of massively parallel sequencing technologies for resolving major radiations in ray-finned 

14 fishes unrealized. Here, we provide a genomic perspective on longstanding questions regard- 

15 ing the diversification of major groups of ray-finned fishes through targeted enrichment of 

16 ultraconserved nuclear DNA elements (UCEs) and their flanking sequence. Our workflow 
l? efficiently and economically generates data sets that are orders of magnitude larger than 
is those produced by traditional approaches and is well-suited to working with museum speci- 

19 mens. Analysis of the UCE data set recovers a well-supported phylogeny at both shallow and 

20 deep time-scales that supports a monophyletic relationship between Amia and Lepisosteus 

21 (Holostei) and reveals elopomorphs and then osteoglossomorphs to be the earliest diverging 

22 teleost lineages. Divergence time estimation based upon 14 fossil calibrations reveals that 

23 crown teleosts appeared 270 Ma at the end of the Permian and that elopomorphs, osteoglos- 

24 somorphs, ostarioclupeomorphs, and euteleosts diverged from one another by 205 Ma during 

25 the Triassic. Our approach additionally reveals that sequence capture of UCE regions and 

26 their flanking sequence offers enormous potential for resolving phylogenetic relationships 

27 within ray-finned fishes. 
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28 Introduction 



29 The ray-finned fishes (Actinopterygii) constitute the dominant radiation of vertebrates on 

30 the planet including more than 32,000 species and equaling or exceeding richness estimates 

31 for the combined total of birds, mammals, and reptiles. Despite a long history of sys- 

32 tematic study, resolution of phylogenetic relationships within this vast radiation remains 

33 elusive. Studies based upon traditional morphological and single-gene, PCR-based molecu- 

34 lar approaches have succeeded in delineating several major lineages of ray-finned fishes, but 

35 conflict over how these lineages are related to one another remains. For example, the earliest 

36 morphological studies of ray-finned fishes unite gar (Lepisosteus) with the bowfin (Amia) 

37 in the clade Holostei [1,2] though this clade is not recovered in some later analyses [3,4]. 

38 The early branching of teleost lineages has also been historically contentious. Systematists 

39 agree on the four earliest-diverging lineages: the osteoglossomorphs (bony-tongues; arawanas, 

40 elephant fishes, and allies), the elopomorphs (tarpons, bonefishes, and eels), the ostarioclu- 

41 peomorphs (anchovies and herrings, minnows, characins, catfishes, and electric eels), and the 

42 euteleosts (salmons, pikes, lizardfishes, and perch-like fishes). However, there is disagreement 

43 over both the relationships among these groups and the basal divergences within euteleosts. 

44 Recent morphological and molecular studies support a sister-group relationship between os- 

45 tarioclupeomorphs and euteleosts [5-7], but beyond this there is little agreement regarding 

46 the relationship among these ancient teleost lineages. Morphological analyses alternatively 

47 place the osteoglossomorphs [8] or the elopomorphs [5,9-11] as the sister group to all other 

48 teleosts and the remaining lineages sister to the ostarioclupeomorph/euteleost clade. Some 

49 molecular analyses place elopomorphs and osteoglossomorphs as the sister group to remain- 
so ing teleosts [12, 13] while others recover a basal divergence between osteoglossomorphs and 

51 other teleosts [7, 14] . 

52 Recently, Near et al. [14] used wide-spread taxonomic sampling, in conjunction with se- 

53 quence collected from nine commonly used nuclear genes, to provide a more comprehensive 

54 phylogenetic hypothesis of relationships among fishes. Their results supported the mono- 
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55 phyly of the Holostei, suggested that the elopomorphs formed the earliest diverging teleost 

56 lineage [15], and provided a new timescale for the divergence of ray-finned fishes. Although 

57 promising, these new insights into the radiation of actinopteryigians relied upon a relatively 

58 modest number of genomic markers, and the stability and timing of these relationship en- 

59 coded throughout the genomes of the target groups remain largely untested. One exception 
eo to this statement includes a recent study by Zou et al. [16] which used transcriptome se- 
ei quences to examine basal divergences within euteleosts. However, the Zou et al. [16] study 

62 did not include several anciently diverging lineages (e.g. Amia, osteoglossomorphs) informing 

63 questions about the early evolution of major groups of ray-finned fishes. 

64 Phylogenomics and next-generation sequencing technologies offer enormous promise for 

65 resolving relationships within actinopterygians and other major sections of the Tree of Life, 
ee However, revolutions within genomics and informatics have had a surprisingly modest effect 
e? on data collection practices within the phylogenetics community: most studies of non-model 

68 organisms continue to rely upon direct sequencing of a moderate number of loci, and work- 

69 flows that do take advantage of massively parallel sequencing platforms remain bottlenecked 

70 by cross-species amplification of phylogenetically informative loci. Several alternatives to 

71 traditional phylogenetic workflows exist that help to overcome the inefficiencies of gene- 

72 based sequencing. One class of these methods is exemplified by the recent work of Zou et 

73 al. [16], who used a combination of de novo transcriptome sequencing, existing transcript 

74 data, and computational methods to identify 274 orthologous groups from which they in- 

75 ferred the phylogeny of the Actinopterygii. The benefits of their approach include the use of 

76 existing, transcript-related data sets (ESTs in GenBank); reasonably well-established data 

77 generation methods; and the collection of data from hundreds of loci across the genomes of 

78 the focal taxa. Limitations of this approach include reliance on sampling fresh or properly 

79 preserved tissues (generally precluding the use of thousands of existing museum samples), 
so dependence of the approach on expression patterns of the tissue sampled, and collection of 
si data from fewer genomic locations than alternative methodologies. 
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82 A second class of phylogenomic methods involves sequence capture of nuclear regions 

83 flanking and including ultraconserved elements (UCEs) [17]. Rather than sequencing ex- 

84 pressed portions of the genome, the UCE-based approach involves enriching organismal 

85 DNA libraries for hundreds to thousands of UCEs and their flanking regions; sequencing 

86 these libraries using massively parallel sequencing; and assembling, aligning, and analyzing 
s? the resulting data using informatic tools. This approach has been successfully used in mam- 
ss mals [18], birds [17,19], and reptiles [20] to generate phylogenomic data sets that contain at 

89 least one order of magnitude more characters than those generated using PCR and to resolve 

90 historically contentious sections of the Tree of Life [18,20]. The UCE approach differs from 

91 transcript-based phylogenomic studies [16] because data collection is independent of expres- 

92 sion pattern, researchers can prepare and enrich libraries from existing tissue collections, and 

93 UCE loci may be better conserved and more numerous across distantly related taxa [18]. 

94 Here, we apply the UCE approach to ray-finned fishes by developing a novel set of 

95 sequence capture probes targeting almost 500 UCE regions in ray-finned fishes. We use the 

96 UCE data to provide the first phylogenomic perspective based upon widespread sampling of 

97 hundreds of markers across the genome on long-standing controversies regarding relationships 

98 at the base of the ray-finned fish Tree of Life. These include whether Lepisosteus and Amia 

99 form a monophyletic group (the Holostei [1,2,21]) and how the major lineages of teleosts, 

100 which constitute >99% of ray-finned fishes, are related to one another [5-7,9-11,22,23]. We 

101 also use 14 fossil calibrations to provide the first time-scale for ray-finned fishes based upon 

102 UCE regions and their flanking sequence. Our results reveal that sequence capture of UCE 

103 regions can efficiently and economically generate massive data sets with strong resolving 

104 power at both deep and shallow phylogenetic scales within fishes. 
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105 Materials and Methods 



we Identification of UCE regions 

107 To identify ultraconserved elements (UCEs) in fishes, we used genome-to-genome alignments 

los of stickleback (Gasterosteus aculeatus) to medaka (Oryzias latipes) to locate nuclear DNA 

109 regions of 100% conservation greater than 80 bp in length. To enable efficient capture-probe 

no design, we buffered these regions to 180 bp (where needed) by including equal amounts of 

in medaka sequence 5' and 3' to each UCE. We aligned or re- aligned these buffered regions 

n2 to the genome-enabled fishes (zebrafish, Danio rerio, stickleback, medaka, and two species 

n3 of puffers, Tetraodon nigroviridis and Takifugu rubripes) using LASTZ [24], keeping only 

lu non-duplicate matches of > 120 bp and > 80% sequence identity across all species in the 

115 set. Based on the intersection of UCE loci across all fishes that were greater than 10 Kbp 

lie apart, we designed a pilot set of 120 bp sequence capture probes for each of the UCEs 

H7 present among all members of the set by tiling probes at 4X density. We had these probes 

us commercially synthesized into a custom SureSelect target enrichment kit (Agilent, Inc.). We 

119 used a higher than normal [25] tiling density to help ameliorate potential sequence differences 

120 among species introduced by buffering shorter UCEs to 180 bp. 

121 Library preparation, UCE enrichment, sequencing, and assembly 

122 We prepared DNA libraries from 18 fish species, including representatives of five acantho- 

123 morph orders and two families of perciforms (Table 2), by slightly modifying the Nextera 

124 library preparation protocol for solution-based target enrichment [17] and increasing the 

125 number of PCR cycles following the tagmentation reaction to 20. Following library prepa- 

126 ration, we substituted a blocking mix of 500 /xM (each) oligos composed of the forward and 

127 reverse complements of the Nextera adapters for the Agilent-provided adapter blocking mix 

128 (Block #3). We incubated species-specific libraries with synthetic RNA probes from the 

129 SureSelect kit for 24 h at 65°C. We followed the standard SureSelect protocol to enrich DNA 
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130 libraries following hybridization; we eluted clean, enriched DNA in 30 fiL of nuclease free 

131 water; and we used 15 yuL of enriched template in a 50 fiL PCR reaction of 20 cycles com- 

132 bining forward, reverse, and indexing primers with Nextera polymerase to add a custom set 

133 of 24 indexing adapters [26]. PCR clean-up was completed using Agencourt AMPure XP. 

134 We quantified enriched, indexed libraries using qPCR (Kapa Biosystems), and we prepared 

135 two library pools containing 10 libraries at equimolar ratios prior to sequencing. 

we We sequenced each pool of enriched DNA using two lanes of a single-end 100 bp Illu- 

137 mina Genome Analyser (GAIIx) run. After sequencing, we trimmed adapter contamination, 

138 low quality bases, and sequences containing ambiguous base calls using a pipeline we con- 

139 structed (https://github.com/faircloth-lab/illumiprocessor). We assembled reads, 
wo on a species-by-species basis, into contigs using Velvet [27] and VelvetOptimiser. Follow- 

141 ing assembly, we used a software package (https : //github . com/f aircloth-lab/phyluce) 

142 containing a custom Python program (match_contigs_to_probes.py) integrating LASTZ [24] 

143 to align species-specific contigs to the set of probes/UCEs we used for enrichment while 

144 removing reciprocal and non- reciprocal duplicate hits from the data set. During matching, 

145 this program creates a relational database of matches to UCE loci by taxon. This pro- 

146 gram also has the ability to include UCE loci drawn from existing genome sequences, for 

147 the primary purpose of including available data from genome-enabled taxa as outgroups 

148 or to extend taxonomic sampling. We used this feature to include UCE loci we identified 

149 in the genome sequences of Gasterosteus aculeatus, Haplochromis burtoni, Neolamprologus 

150 brichardi, Oreochromis niloticus, Oryzias latipes, Pundamilia nyererei, Takifugu rubripes, 

151 Tetraodon nigroviridis, Gadus morhua, and Lepisosteus oculatus. After generating the rela- 

152 tional database of matches to enriched sequences and genome-enabled taxa, we used addi- 

153 tional components of PHYLUCE (get_match_counts.py) to query the database and generate 

154 fasta files for the UCE loci we identified across all taxa. Then, we used a custom python 

155 program (seqcap_align_2.py) to align contigs with MAFFT [28] and trim contigs representing 

156 UCEs, in parallel, across the selected taxa prior to phylogenetic analysis [17]. 
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157 Phylogenetic Analyses 

158 The large number of UCE loci we collected create a vast potential space for partitioning data 

159 that makes a traditional evaluation of alternative partitioning strategies computationally in- 
leo tractable. As a result, we modeled nucleotide substitutions across the concatenated data 
lei set using two approaches. For Bayesian analysis, we used a custom script (runjnraic.py) 

162 wrapping a modified MrAIC 1.4.4 [29] to find the best-fitting model for each UCE locus, we 

163 grouped loci having similar substitution models (selected by AICc) into the same partition, 

164 and we assigned the partition specific substitution model to all loci concatenated within each 

165 partition. For maximum likelihood analyses, we maintained the partitions identified in the 
lee Bayesian analysis and we modeled each partition using the GTR+CAT approximation. We 

167 performed Bayesian analysis of the concatenated data set using MrBayes 3.1 [30] and two 

168 independent runs (4 chains each) of 5,000,000 iterations each, sampling trees every 500 iter- 

169 ations, to yield a total of 10,000 trees. We sampled the last 5,000 trees after checking results 

170 for convergence by visualizing the log of posterior probability within and between the inde- 

171 pendent runs for each analysis, ensuring the average standard deviation of split frequencies 

172 was < 0.001, and ensuring the potential scale reduction factor for estimated parameters was 

173 approximately 1.0. We performed maximum likelihood analysis of the concatenated data in 

174 RAxML [31] using the rapid bootstrapping algorithm and 500 bootstrap replicates. 

175 Gene tree-species tree methods enjoy some advantages over the analysis of concatenated 

176 data sets under certain conditions [32-34] but may also be sensitive to missing data [35] and 

177 to the resolution of individual gene trees [36]. To minimize the number of unresolved gene 

178 tree topologies and maximize the number of topologies that overlapped in sampling the base 

179 of the actinopterygian tree, we selected a subset of the UCE contigs containing complete 

180 data for Polypterus and Acipenser and loci > 50 bp, and we used this subset to estimate a 
lei species tree with CloudForest (https://github.com/ngcrawford/CloudForest), a parallel 

182 implementation of a workflow combining substitution model selection (identical to MrAIC 

183 1.4.4 [29]) and genetree estimation using PhyML [37]. We estimated the species tree by 
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184 summarizing gene trees using STAR [38-40]. To assess confidence in the resulting species 

185 tree, we used CloudForest to generate 1000, multi-locus, non-parametric bootstrap replicates 
we by resampling nucleotides within loci as well as resampling loci within the data set [41] , we 
187 summarized bootstrap replicates using STAR, and we reconciled bootstrap replicates with 
las the species tree using RAxML. 

189 Divergence Time Estimation 

190 We used a set of 14 calibration points (Appendix: Fossil calibrations) from previous timetree 

191 studies [42, 43] to date several key splits on the tree ranging in age from the Givetian (392 

192 Ma) to the early Rupelian (32 Ma). We used the RAxML topology (Fig. 1) plus fossil- 

193 derived minimum and maximum age constraints to infer divergence times in BEAST vl.72 

194 [44] assuming a GTR model of sequence evolution with gamma-distributed rate variation. 

195 Preliminary analyses with the full UCE data (all 491 loci) showed poor mixing under a wide 
we range of fossil constraint parameterizations, and all models incorporating the full UCE data 
197 set failed to reach convergence after 300 million generations. To estimate divergence times 
we for these taxa from a large data set, we randomly sampled 50 UCE loci, with replacement, to 

199 create five different matrices for BEAST analyses. We analyzed the resulting matrices using 

200 independent runs of BEAST. MCMC runs of these data sets mixed much better and showed 

201 strong evidence of convergence, having ESS values > 200 for nearly all parameters after 50 

202 million generations. However, the parameters for the mean and variance of the lognormally 

203 distributed rates mixed much more slowly with autocorrelation times of > 5* 10 5 generations 

204 and yielding ESS of ~ 50 - 100 across the replicates. 
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205 Results and Discussion 



206 Probe design, UCE enrichment, and sequencing 

207 We located 500 UCEs shared among all actinopterygian fishes, and we designed a set of 

208 2,000 capture probes targeting each of these loci (4X coverage). Following enrichment and 

209 sequencing, we obtained an average of 2,819,047 reads per species, which we assembled into an 

210 average of 665 contigs having an average length of 457 bp (Table 2). After removing contigs 

211 that matched no UCEs and UCE loci that matched multiple contigs, we enriched an average 

212 of 332 (50%) unique contigs matching UCE loci from each species. Average sequencing depth 

213 across UCE loci was 498X. We integrated extant genomic data from several fish species to 

214 this group of unique UCE contigs, and we constructed 491 alignments (average length: 305 

215 bp, 95 CI: ±16.0) comprised of 149,246 characters. Each alignment contained an average of 

216 21 target taxa (95 CI ± 0.4) after data trimming, and we used this incomplete data matrix 

217 for subsequent analyses with RAxML and MrBayes. After removing loci having missing data 

218 for Polypterus and Acipenser, we input 136 alignments (41,731 characters; average length: 

219 307 bp, 95 CI: ± 27.7) to CLOUDFOREST for model selection and subsequent species tree 

220 estimation using STAR. 

221 A phylogenomic perspective on the basal radiation of ray-finned 

222 fishes 

223 Maximum likelihood analysis produced a single, completely resolved topology wherein all but 

224 two nodes received high (> 0.99) bootstrap proportions and Bayesian posterior probabilities 

225 (Fig. 1). This topology provides new insight into several long-standing questions concerning 

226 the evolution of ray-finned fishes. Our analysis strongly supports the monophyly of the 

227 Holostei (Amia + Lepisoteus) . This clade is historically controversial because morphological 

228 studies alternatively support [1,2,21] and refute [3,4] the monophyly of this group, while 

229 recent molecular studies generally recover the relationship [15,45,46]. Additionally, our 
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230 analyses do not support prior findings of an "ancient fish clade" including the Holostei + 

231 Aciper serif ormes as the sister group to the teleosts [46,47]. Rather, our results strongly 

232 suggest a traditional relationship in which these lineages form successive sister groups to the 

233 teleosts. 

234 Our phylogenomic data provide strong evidence for the placement of elopomorphs as 

235 the sister group to all other teleosts and osteoglossomorphs and ostarioclupeomorphs as 

236 successive sister lineages to the euteleosts (Fig. 1). Our maximum likelihood topology is 

237 strongly incongruent with mitogenomic studies [7, 14] but consistent with both a recent 

238 analysis of multiple nuclear genes [15] and some of the earliest morphological analyses of 

239 the group [5, 9-11]. Within euteleosts, our results are congruent with recent molecular 

240 studies [6,15, 16] in placing esociforms as the sister to salmoniforms rather than any neoteleost 

241 lineages. 

242 Within acanthomorphs, the largest clade of euteleosts, UCEs recover several intriguing 

243 clades that agree with results from recent molecular phylogenetic studies. These include 

244 the African cichlids + medaka (Clade CI, Fig. 1), corresponding to an expanded clade of 

245 atherinomorphs suggested by recent studies [16,48,49]; a clade of gasterosteiforms (stickle- 

246 back) and scorpaeniforms (Taenionotus) that is congruent with recent molecular and mor- 
24? phological studies [16,50,51]; and a clade including surgeonfish, frogfishes, and pufferfishes 

248 (acanthuroids, lophiiforms, and tetraodontiforms) corresponding to acanthomorph clade "N" 

249 of Dettai and Lecointre [15,52]. Together with the estimated topology, our timetree suggests 

250 that UCEs provide sufficient phylogenetic signal to resolve divergences within fishes that are 

251 as young as 5 Ma and as old as 400 Ma. 

252 The STAR topology was less resolved than topologies based upon analysis of the con- 

253 catenated data set (Fig. SI) but recovered largely congruent relationships including a mono- 

254 phyletic Holostei as the sister to other actinopterygians, monophyly of elopomorphs, os- 

255 teoglossomorphs, ostarioclupeomorphs, and euteleosts, and a successive sister group rela- 

256 tionship between ostarioclupeomorphs, Salvelinus + Umbra, and all remaining euteleosts. 



11 



257 The species tree switched the position of the Gadiformes, represented by cod (Gadus) and 

258 Myctophiformes, represented by Diaphus. This resolution is not congruent with results from 

259 Near et al, [15] but has been suggested before in other molecular studies [6,53]. 

260 A phylogenomic timescale for the radiation of ancient fish lineages 

261 A general expectation during Bayesian analysis is that runs will converge more quickly on the 

262 true posterior distribution with increasing data. However, in our divergence time analysis, 

263 we were only able to obtain acceptable mixing in BEAST when analyzing subsamples of 

264 the data set; analyses of all 491 loci did not converge despite a diversity of partitioning and 

265 calibration strategies combined with long MCMC analyses. This was surprising given the 

266 relatively modest number of taxa included in this study and because we fully constrained the 

267 topology. Discordant fossil calibrations [54] could potentially underlie poor mixing behavior. 

268 However, alternative fossil calibration strategies that enforced subsets of the 14 calibrations 

269 revealed that even a single constraint on the age of the root node resulted in analyses with 

270 poor mixing behavior. We observed additional evidence that fossil conflict was not the 

271 cause of poor mixing because subsampled analyses including all 14 constraints converged 

272 and produced nearly identical credible intervals for all nodes (Fig. S2). One possible barrier 

273 to good mixing during analyses of the full data set in BEAST may be that the proposal 

274 mechanisms and/or tuning parameters for branch lengths are not suited for alignments that 

275 contain both very quickly (UCE flanks) and very slowly (UCE core regions) evolving sites. 

276 With a large number of sites, even a relatively small change in branch length could result in 

277 a large difference in likelihood, resulting in a high rejection rate of branch length proposals. 

278 This explanation is consistent with the extremely low acceptance rates we observed in our 

279 analyses ( < 5% for branch rate and node height proposals). The recovery of similar age 

280 estimates across UCE subsamples (Fig. S2) suggests that the estimated time scale for ray- 

281 finned fish diversification adequately reflects the signal for divergences contained within the 

282 full UCE data set. 
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283 Our divergence time analyses suggest that the extant radiation of ray-finned fishes began 

284 rs 420 Ma near the end of the Silurian. Divergences amongst the actinopterians, neoptery- 

285 gians, and holosteans (Fig. S2, nodes 1-4) are very ancient and span the late Devonian to 

286 the early Permian. Extant teleosts (Fig. S2, node 5) can trace their origin back to ~ 270 

287 Ma during the late Permian with primary divergences amongst the elopomorphs, osteoglos- 

288 somorphs, ostarioclupeomorphs, and euteleosts occurring before the end of the Triassic. The 

289 credible intervals surrounding the age of teleosts are wide, however, and teleosts may have 

290 originated as early as 309 Ma during the upper Carboniferous to as late as 226 Ma during 

291 the mid- Triassic. We recover a basal split between Salvelinus + Umbra, representing two 

292 protancanthopterygian lineages, and the rest of the euteleosts, at 170 Ma during the Juras- 

293 sic (Fig. S2, node 9). Extant acanthomorphs, which comprise the bulk of teleost diversity 

294 including over 16,000 species and 300 families, trace their origins to the Cretaceous, 124 

295 Ma (Fig. S2, node 14). We report ages for all splits in Table 1. 

296 Our analysis represents the first timescale for ray-finned fishes derived from widespread 

297 sampling of the nuclear genome. Previous divergence time analyses based upon fossils, mi- 

298 tochondrial DNA, and nuclear DNA have produced conflicting timescales for ray-finned fish 

299 diversification. For example, the earliest known fossil teleosts are elopomorphs and ostario- 

300 physans, dating to the Late Jurassic [42,43]. In contrast mitogenomic data consistently sug- 

301 gest Paleozoic divergences, usually into the early Permian or late Carboniferous (e.g., [55]). 

302 Some divergence time studies based upon nuclear genes have posited Late Triassic to Early 

303 Jurassic divergences (173-214 Ma [42,43]) for the teleosts, while more recent studies that 

304 include additional nuclear loci and an improved set of fossil calibrations push the age of 

305 crown teleosts to 307 Ma at the end of the Carboniferous [15]. Our UCE-derived date for 

306 the divergence of teleosts (~ 270 Ma) is much older than an earlier estimate based upon 

307 RAG1 sequence [43] and largely overlaps with the age estimate derived by Near et al [15]. 

308 Other splits within our UCE timescale generally correspond to those derived by Near et 

309 al [15] (Fig. S3). One major exception to this observation is divergence of the euteleosts 
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310 which UCE data estimate to be 171 Ma (95% credible interval, 152-194 Ma) while tra- 

311 ditional nuclear loci estimate an age of ~ 229 Ma (95% credible interval 220-259 Ma) [15]. 

312 One possible cause for this difference is age estimates is that we have not sampled some of 

313 the earliest diverging eutelost lineages in our study including Lepidogalaxias, resulting in an 

314 underestimate of the true euteleost crown age. Our molecular timescale provides additional 

315 evidence for a relatively ancient origin of teleosts and further highlights the apparent gap in 

316 the fossil record between stem and crown members of this clade [15,42]. 

317 Conclusions 

318 Sequence capture of regions anchored by UCEs offers a powerful and efficient means of gener- 

319 ating massive genomic data sets capable of resolving phylogenetic relationships at both deep 

320 and shallow scales in non-model organisms. Our UCE-based approach offers several advan- 

321 tages over previous studies that should contribute to the reliability of our topology. These 

322 benefits include efficient sampling of sequence data across individual genomes and among 

323 divergent taxa, collection of data from an order of magnitude more loci than studies based 

324 upon traditionally-used genetic markers and almost twice as many loci as transcriptome- 

325 based genomic studies [16], validity of the UCE probe set across bony fishes spanning 350 

326 Ma of evolutionary history, and utility of the UCE enrichment approach with tissues col- 

327 lected from museum specimens. Additionally, these data illustrate that biologists can use 

328 UCE-based genetic markers to reconstruct the phylogeny of taxa other than amniotes, sup- 

329 porting the observation that UCE-based markers are a universal source of phylogenetically 

330 informative characters [17,18]. 
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bib Description of timetree calibration points 

516 We assigned fossil calibrations to 14 nodes of the maximum likelihood phylogeny (Fig. 2) 

si? to enable divergence time estimation using BEAST vl.72. Below, we justify the bounds for 

sis each calibration point. 



519 MRCA of Actinopterygii 

520 Fig. 2, node 1; uniform calibration = 392-472 Ma 

521 We used the stegotrachelids fossils from the Givetian/Eifelian boundary (^ 392 Ma) 

522 used by [42] to date the minimum age of the root and assumed an upper bound of 

523 472 Ma based upon the minimum age for the split between acanthodians and all other 

524 bony fishes [56]. 

525 MRCA of Actinopteri 

526 Fig. 2, node 2, exponential calibration lower bound = 345 Ma, upper 95% 

527 = 392 Ma 

528 The oldest fossil belonging to this clade is Cosmoptycius from the Tournasian (Car- 

529 boniferous, 359-345 Ma) [42]. We used the stegotrachelids fossils from the Give- 

530 tian/Eifelian boundary (« 392 Ma) used by [42] to date the crown actinopterygians to 

531 set the upper bound. 

532 MRCA of Holostei 

533 Fig. 2, node 3; exponential calibration lower bound = 284 Ma, upper 95% 

534 = 345 Ma 

535 The oldest fossil assigned to this clade is the neopterygian Brachydegma caelatum, 

536 from the Artkinsian (early Permian, 284 Ma) [42]. We used Cosmoptychius from the 

537 Tournasian (Carboniferous, 359-345 Ma) to set the upper bound [42]. 

538 MRCA of Teleostei (Fig. 2, node 5; unconstrained) 

539 The oldest crown teleost is Anaethalion, from the late Kimmeridgian (Jurassic, 152 Ma) 
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540 [43]. However, due to the appearance of fossils from several major teleost groups (e.g., 

541 the ostarioclupeomorph Tischlingerichthys viohli, the euteleost Leptolepides spratti- 

542 formis) in deposits of the same age it seems likely that teleosts arose considerably 

543 earlier. We left this node unconstrained in our analysis to obtain an estimate based 

544 upon other actionpterygian calibrations and the UCE sequence data. 

545 MRCA of Osteoglossomorpha 

546 Fig. 2, node 13; lower bound = 112 Ma, upper 95% = 225 Ma 

547 The oldest taxon assigned to our osteoglossomorph clade (which only includes a sub- 

548 set of the major osteoglossomorph lineages) is the arapamid Laeliichthys ancestralis 

549 from the Aptian Areado Formation of Brazil (Cretaceous, 112 Ma) [57]. We used the 

550 Hiodontidae Yanbiania wangqingica from the Barremian (Cretaceous, 130 Ma) to set 

551 the upper boundary [43]. 

552 MRCA of Ostarioclupeomorpha 

553 Fig. 2, node 8; lower bound = 149 Ma, upper 95% = 225 Ma 

554 The oldest taxon assigned to this clade is Tischlingerichthys viohli, from the upper 

555 Tithonian (Jurassic, 149 Ma) [43]. We used the stem teleost Pholidophoretes salvus 

556 (Pholidophoridae), from the early Carnian/ Julian (Triassic, 228-225 Ma) to set the 

557 upper bound [43]. 

558 MRCA of Euteleostei 

559 Fig. 2, node 9, lower bound = 152 Ma, upper 95% = 225 Ma 

560 The oldest taxon assigned to this clade is Leptolepides sprattiformis (Orthogoniklei- 

561 thridae) from the late Kimmeridgian (Jurassic, 152 Ma) [43]. We used the stem teleost 

562 Pholidophoretes salvus (Pholidophoridae), from the early Carnian/ Julian (Triassic, 

563 228-225 Ma) to set the upper bound [43]. 

564 MRCA of Elopomorpha 

565 Fig. 2, node 10, lower bound = 135 Ma, upper 95% = 225 Ma 
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566 The oldest taxon assigned to this clade is the albulid Albuloideorum ventralis, from 

se? the early Hauterivian (Jurassic/Cretaceous border, 135 Ma) [43]. The oldest stem 

see elopomorph is Anaethalion, from the late Kimmeridgian (Jurassic, 152 Ma) [43]. 

569 Characiformes vs Cypriniformes 

570 Fig. 2, node 16, lower bound = 100 Ma, upper 95% = 149 Ma 

571 The oldest taxon assigned to this clade is the characiform Santanichthys diasii from 

572 the Albian (Cretaceous, 112-100 Ma) [43]. We used Tischlingerichthys viohli, from the 

573 upper Tithonian (Jurassic, 149 Ma) to set the upper bound [43]. Our prior assumed 

574 100 My as the minimum age, and 149 Ma for the upper bound. 

575 Salmoniformes vs Esociformes 

576 Fig. 2, node 12, lower bound = 125 Ma, upper 95% = 152 Ma 

577 The oldest taxon assigned to this clade is Helgolandichthys schmidi, from the early 

578 Aptian (early Cretaceous, 125 Ma) [43]. We used Leptolepides sprattiformis (Or- 

579 thogonikleithridae) from the late Kimmeridgian (Jurassic, 152 Ma) to set the upper 

580 bound [43]. 

581 MRCA of Ctenosquamata 

582 Fig. 2, node 11, lower bound = 122 Ma, upper 95% = 152 Ma 

583 The otoliths assigned to Acanthomorphorum forcallensis from the Aptian (Early Cre- 

584 taceous, 124-122 Ma) represents the oldest fossils assigned to this clade [43]. We used 

585 Leptolepides sprattiformis (Orthogonikleithridae) from the late Kimmeridgian (Juras- 
sse sic, 152 Ma) to set the upper bound [43]. 

587 MRCA of Acanthopterygii 

sss Fig. 2, node 15, lower bound = 99 Ma, upper 95% = 122 Ma 

589 The oldest taxa assigned to this clade are various Beryciformes (e.g., Hoplopteryx, 

590 Trachichthyoides) from the Cenomanian (Late Cretaceous, 99 Ma) [43]. We used the 
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591 fossil otoliths assigned to " Acanthomorphorum" forcallensis from the Aptian (Early 

592 Cretaceous, 124-122 Ma) to set the upper bound [43]. 

593 Gasterosteus vs Taenianotus 

594 Fig. 2, node 21, lower bound = 85 Ma, upper 95% = 122 Ma We used the 

595 oldest gasterosteiform, Gasterorhamphus zuppichinii from the Santonian (84-85 Ma) 

596 to date the minimum age of this split [57]. We used the fossil otoliths assigned to 

597 "Acanthomorphorum" forcallensis from the Aptian (Early Cretaceous, 124-122 Ma) to 

598 set the upper bound [43]. 

599 Lophiiforms vs tetraodontiforms 

600 Fig. 2, node 20, lower bound = 85 Ma, upper 95% = 122 Ma 

601 The oldest taxon assigned to this clade is the stem tetraodontiform Cretatriacanthus 

602 guidottii from the Santonian of Nardo (Italy). We chose this taxon to date the min- 

603 imum age rather than other, older, stem tetraodontiformes because preliminary re- 

604 examination of the relationships of extant and fossil tetraodontiforms (Santini, unpub- 

605 lished) casts doubt on their phylogenetic affinities. We used the fossil otoliths assigned 

606 to " Acanthomorphorum" forcallensis from the Aptian (Early Cretaceous, 124-122 Ma) 
60? to set the upper bound [43]. 

608 MRCA of Tetraodontidae 

609 Fig. 2, node 23, lower bound = 32 Ma, upper 95% = 50 Ma 

610 The oldest taxon assigned to this clade is Archaeotetraodon winterbottomi from the 
en Oligocene of Caucasus (32-35 Ma) [43]. We used the stem tetraodontid Eotetraodon 
6i2 pygmaeus from the Ypresian (middle Eocene, 50 Ma) to set the upper bound [58]. 
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Figure Legends 



■ Polypterus 




Figure 1. Maximum likelihood phylogram of ray-finned fish relationships based 
upon UCE sequences. All nodes except for two (indicated by arrows) supported by 
bootstrap proportions and Bayesian posterior probabilities > 0.99. Our analysis supports a 
monophyletic Holostei and reveals the elopomorphs to be the earliest diverging lineage of 
teleosts. CI, C2, and C3 indicate clades within acanthomorphs consistent with other recent 
molecular studies (see Discussion). 
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Figure 2. UCE-based timescale of ray-finned fish diversification. Arrows indicate 
fossil-calibrated nodes (described in Appendix 1). Numbers refer to splits in Fig. 1. 



29 



Polypterus 
Acipenser 



100 



100 



100 



f I — 

I— I 



100 



100 _ 



100 



100 



100 



Lepisosteus 
Amia 

Strophidon 

Megalops 

Osteoglossum 

Pantodon 

Anchoa 

Danio 

Astyanax 

Salvelinus 

Umbra 

Gadus 

Diaphus 

Lampris 

Antennarius 

Takifugu 

Tetraodon 

Oryzias 

Taenianotus 

Gasterosteus 

Acanthurus 

Pundamilia 

Neolamprologus 

Haplochromis 

Oreochromis 



Figure SI. Species tree based upon STAR analysis. Topology based upon analysis 
of all loci > 50 base pairs that contained both Polypterus and Acipenser (N = 136). Node 
values indicate bootstrap proportion based upon 1000 replicates. We collapsed nodes 
having < 50% bootstrap support. 




Node 



Figure S2. Mean and 95% credible interval of divergence time estimates for subsamples of 
the UCE data set. Node numbers refer to labels in Fig. 2. Arrows indicate 
fossil-constrained nodes. For each node, we show the mean and 95% credible interval of five 
divergence time analyses. 
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Figure S3. Comparison of ages derived from UCE data and nine nuclear genes [15]. UCE 
ages are in green. Circles indicate mean age and error bars indicate 95% credible interval. 



32 



Tables 



Table 1. Divergence times. Node numbers refer to Fig. 2. We averaged mean and 95% 
credible interval estimates (in Ma) across five subsampled divergence time analyses (see 
text). 



Node 


Clade 


Split 


Age 


Min 


Max 


1 


Actinopterygii 


Polypterus vs Takifugu 


420.0 


392.0 


461.0 


2 


Actinopteri 


Acipenser vs Takifugu 


366.5 


345.0 


398.8 


3 


Neopterygii 


Lepisosteus vs Takifugu 


327.1 


295.5 


360.8 


4 


Holostei 


Lepisosteus vs Amia 


295.9 


284.0 


317.1 


5 


Teleostei 


Megalops vs Takifugu 


268.8 


226.8 


309.4 


6 


Osteoglossocephala 


Pantodon vs Takifugu 


251.7 


211.1 


291.7 


7 


Clupeocephala 


Anchoa vs Takifugu 


210.9 


177.1 


245.3 


8 


Ostarioclupeomorpha 


Anchoa vs Astyanax 


169.5 


149.0 


197.8 


9 


Euteleostei 


Umbra vs Takifugu 


170.6 


152.0 


193.7 


10 


Elopomorpha 


Megalops vs Strophidon 


160.9 


135.0 


206.6 


11 


Ctenosquamata 


Diaphus vs Takifugu 


132.1 


122.0 


145.8 


12 


protacanthopterygii 


Umbra vs Salvelinus 


130.9 


125.0 


141.8 


13 


Osteoglossiformes 


Pantodon vs Osteoglossum 


128.9 


112.0 


161.5 


14 


Acanthomorphs 


Gadus vs Takifugu 


125.4 


112.2 


140.6 


15 




Lampris vs Takifugu 


116.4 


104.0 


129.5 


16 


Ostariophysans 


Danio vs Astyanax 


111.4 


100.0 


131.9 


17 




Oryzias vs Takifugu 


104.1 


94.5 


114.6 


18 




Gasterosteus vs Takifugu 


99.1 


91.0 


108.2 


19 




Acanthurus vs Takifugu 


94.1 


87.4 


102.0 


20 




Gasterosteus vs Taenianotus 


89.7 


85.0 


97.0 


21 




Antennarius vs Takifugu 


87.9 


85.0 


93.1 


22 




Oryzias vs Haplochromis 


76.7 


42.1 


103.3 


23 


Tetraodontidae 


Tetraodon vs Takifugu 


37.5 


32.0 


48.0 


24 


Cichlidae 


Oreochromis vs Haplochromis 


9.9 


3.2 


19.2 


25 




Neolamprologus vs Haplochromis 


4.7 


1.3 


9.4 


26 




Pundamilia vs Haplochromis 


2.0 


0.3 


4.4 
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